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1.  INTRODUCTION 


Oun  tube  heating  from  multiple  firings  continues  to  be  a  subject  of  concern  to  ordnance  engineers. 
A  number  of  investigators  have  modeled  the  heating.  Some  recent  publications  include:  Artus  and 
Hasenbein  (1989):  Chandra  and  Fisher  (1989a.  1989b):  Talley  (1989a.  1989b):  Rapp  (1990):  Chandra 
(1990):  Qerber  and  Bundy  (1991):  and  Conroy  (1991).  In  these  studies,  the  thermal  conductivity,  density, 
and  specifle  heat  of  the  metal  were  taken  to  be  constants,  and  the  resulting  problems  were  linear. 
However,  the  conductivity  and  specific  heat  are,  in  fact,  functions  of  temperature. 

During  the  heat  transfer  process  from  a  single  gun  firing,  the  temperature  of  the  metal  near  the  bore 
surface  can  rise  more  than  SOO  K.  The  specific  heat  and  thermal  conductivity  will  vary  significantly  over 
such  a  large  temperature  interval.  According  to  Fourier's  heat  flux  principle,  we  expect  an  increase  in  the 
conductivity  with  temperature  to  augment  the  resulting  heat  transfer,  and  hence  further  Increase  the 
temperature  over  a  scenario  in  which  the  conductivity  is  constant.  On  the  other  hand,  an  increase  in 
specific  heat  with  temperature  will  have  the  opposite  effect— it  will  tend  to  diminish  the  subsequent 
temperature  rise.  Thus,  it  is  difficult  to  discern,  a  priori,  what  effect  temperature-dependent  thermal 
properties  will  have  on  gun  tube  heating.  We  propose  to  answer  this  question  by  solving  the  nonlinear 
heat  transfer  problem  with  a  finite  difference  solution,  using  the  Ciank-Nicolson  implicit  scheme  with 
iteration.  Application  of  the  Kirchhoff  transform^on,  described  in  Section  3,  will  render  the  problem 
more  tractable. 

We  choose  to  apply  the  model  to  gun  barrel  heating  in  uink  guns,  specifically,  the  120-mm  M2S6 
cannon  for  the  Ml  A1  tank.  A  tank  gun  is  ammunition-limited  by  the  onboard  stowage  restrictions.  This 
imparts  an  upper  limit  to  the  firing-induced  barrel  temperature  change  for  tank  guns  that  does  not  hold 
for  unlimited- fire  guns,  such  as  small-caliber  chain  guns  or  large-caliber  artillery  guns.  Consequently,  w  e 
do  not  treat  the  broadest  possible  temperature  range  for  all  types  of  guns:  neverthele.ss,  our  findings  cover 
a  substantial  temperature  span.  The  findings  indicate  that  the  Inclusion  of  temperature  dependence  in  the 
modeled  thermal  properties  docs  not  have  a  significant  Impact  on  the  barrel  temperature  predictions  over 
the  range  Investigated. 
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2.  THI3  HEATING  MODEL 


yft  shall  compute  the  gun  barrel  heating  for  multiple  firings  by  an  extension  of  the  model  used  in 
Oeiber  and  Bundy  (1991).  The  following  assumptions  apply  here: 

(1)  Temperature  gradients  in  the  axial  direction  are  neglected  in  comparison  to  those  in  the  radial 
direction. 

(2)  Temperature  is  axisymmetric  in  the  plane  nonnal  to  the  bore  axis.  This  implies  axisynunctric  heat 
input,  as  well  as  the  neglect  of  gravity,  barrel  thickness  variation,  and  other  effects  tliat  would  cause 
azimuthal  dependence. 

(3)  Feedback  of  batrel  heating  to  flow  in  the  gun  bore  is  neglected,  so  that  the  same  bore  temperamre 
and  convective  heat  transfer  coefficient  histories  (for  a  single  round)  furnish  the  input  data  for  every  round 
calculated. 

(4)  Friction  heating  is  neglected. 

(5)  Thermal  expansion  of  the  barrel  is  not  considered  to  have  an  effect  on  the  heat  transfer  process, 

(6)  The  density,  p,  of  the  gun  bairel  metal  is  constant  (>  7,827  kg^i’). 

(7)  The  thermal  conductivity,  k,  and  the  specific  heat,  Cp,  of  the  metal  are  functions  of  the 
temperature,  T,  as  shown  in  Figures  1  and  2,  respectively.  (Figure  1  is  based  on  data  from  the  Mateilal 
Propenies  r^ata  Center  [1973]:  figure  2  is  based  on  data  provided  by  Joseph  Cox,  Benet  Weapons 
Laboratory  [1990].)  Also  Indicated  in  the  figures  are  the  constant  values  used  by  the  authors  for  k  and 
Cp  in  a  previous  study  (Geitcr  and  Bundy  1991):  k  *  38.07  J/(m  s  K)  and  Cp  ■  469.05  JAg  K.  The 
particular  constants  chosen  are  within  the  range  of  values  spanned  in  the  course  of  a  typical  firing 
scenario. 

The  abmpt  changes  in  the  cur\'e  of  Figure  2  are  associated  with  the  phase  transiuon  of  steel  from  a 
body-centered  crystal  lattice  structure  (ferrite)  to  a  face-centered  lattice  structure  (austenite).  Even  though 
we  adjust  the  k  and  Cp  to  account  for  phase  transition  in  the  conduction  eqi.ation,  we  shall  not  (in  this 
report)  further  modify  the  conduction  equation  to  account  foi  the  latent  heat  effects  of  the  phase  transition. 
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3.  FORMULATION  OF  THE  PROBLEM 


3.1  Stflfwnwit  rif  the  Problem.  We  state  our  piublem  in  tenns  of  cylindrical  coordinates:  r,  8,  andz. 
The  radial  coordinate,  r,  is  zero  on  the  axis  of  the  gun  tube  (z-axis),  and  varies  from  R(  to  R,,  the  radii 
of  the  concentric  inner  and  outer  walls,  respectively  (Figure  3).  As  implied  in  Section  2,  the  azimuthal 
angle,  8,  does  not  enter  the  problem.  The  axial  coordinate,  z,  is  taken  to  be  zero  at  the  gun's  breech.  The 
barrel  temperature,  T(r,z,t),  where  t  is  time  measured  from  the  initiation  of  the  first  round,  is  determined 
by  the  following  differential  equsMon  of  heat  conduction  for  a  stationaiy,  homogeneous,  isotropic  solid 
with  no  internal  heat  generation  (Ozisik  1968,  3S3,  or  see  ADDENDUM): 

p  Cp  dT/dt  ■  div  [k  grad  T] . 

Under  the  assumptions  made  in  Section  2,  this  equation  reduces  to 

(1/a)  3T/8t  -  [a*T/ai*  +  (1/r)  ar/ar]  +  (dk/dD  (3T/ar)»/k  ,  (1) 

where  a  ■  k/(p  cp  is  the  theimal  difhisivity,  now  a  ftinction  of  T  as  shown  in  Figure  4. 

The  initial  and  boundary  conditions  do  not  change  in  form  from  those  for  constant  properties.  Let 
T.  designate  the  ambient  temperature  of  the  atmosphere,  assumed  to  be  constant.  The  initial  condition  is 

T(r)  =  T_  ,  t  ■  0.  R,  S  r  i  R,  (z  ■  const) .  (2) 

The  boundary  condition  at  the  inner  wall  is  Newton’s  law  of  cooling: 

-  k  9T^  -  hgCrg-T),  r«R,,t>0  (z»  const) ,  (3) 

where  T,  (t,z)  is  the  cross-sectional  average  temperature  of  the  flow  in  the  bore,  and  h,  (t,z)  is  the 
coefficient  of  heat  transfer  oetween  the  gas-particle  mixture  in  the  bore  and  the  inner  wall  of  the  barrel. 
T,  (u)  and  h,  (t,z)  are  known  flpom  Interior  ballistic  computations  and  thus  constitute  input. 

The  boundary  condition  at  the  outer  wall,  which  includes  both  convective  and  radiative  cooling  (Ozisik 
1968:  Equations  1-28  and  8-1 37c),  is 
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-  k  ar/dr  =  h_  CT  ~  TJ  +  Fo  (T<  -  Tj)  r  =  R,,  t  >  0  (z  =  const),  (4) 

where  >•  h.(z)  Is  the  coeffldent  of  convective  heat  transfer  between  the  barrel  wall  and  the  surrounding 
atmosphere.  F  is  the  radiation  interchange  factor  between  the  barrel  outer  wall  and  the  environment  (in 
our  case,  assumed  *  0.9S),  and  o  is  the  Stefan-Boltzmann  constant  [>  5.669  x  10^  J/(m*  s  K*)]. 

3.2  Kirchhoff  Transfomaation.  Equations  1,  3.  and  4  show  that  the  heat  conduction  problem  is 
nonlinear  in  T.  Introducing  the  Kirchhoff  transformation  (Ozlsik  1968,  353;  Boley  and  Weiner  1960, 
141), 


U<T)-  f^^[k(t)yic,]dt. 


(5) 


will  simplify  the  conduction  equation  and  boundary  condidons.  U  has  the  dimensions  of  temperature. 
Here,  Ti  is  an  arbitrarily  chosen  lower  Integradon  limit,  namely,  a  data  point  in  our  table  for  k; 
T|  m  33.4512  K  and  ki «  kCTt)  ■  12.801  J/(m  s  K).  We  shall  solve  the  problem  posed  in  terms  of  U  and 
then  evaluate  T(U)  from  the  Inverse  of  the  Kirchhoff  transformadon. 

The  variable  coefficients  k(D,  c,(D.  and  a(D  are  available  in  tabular  form.  For  a  particular 
aigiunent,  the  ftincdon  value  is  evaluated  by  3-point  interpoladon.  UCT)  is  obtained  in  tabular  form  by 
numerical  integradon;  it  is  the  solid  curve  in  Figure  5.*  The  problem  now  has  the  form 


(1/a)  au/dt «  3*U/ai^  +  (1/r)  dU/ar  (6a) 

-k,  au/ar  * h,  [T, - T(U)1  r-R,,t>0  (6b) 

-k,  au/ar - h.  (T - TJ  f  FO CT* -  T«*)  r-R,,t>0  (6c) 

U-U.-UCTJ  tSO,  R,SrSR,  (6d) 


*  It  wu  neouiuy  to  txloid  th«  tablu  of  k  md  U  in  Figttru  1  and  S  by  mmpolition  to  higher  temperiturei  (2,000  K)  beeiuie 
of  extremely  high  lurface  tempertturei  etubmd  briefly  in  the  chamber. 
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Figure  5.  Kirehhoff  Transform.  U.  vs.  Temperature. 


The  difAisivity.  a  ■  a[T(U)],  is  a  function  of  U.  shown  in  Hgure  6;  thus.  Equation  6a  is  still  nonlinear. 
However,  the  highly  nonlinear  term  containing  k  «ul  its  derivative  no  longer  appears.  In  order  to  render 
the  flnite-dlffetence  algorithms  for  the  boundary  conditions  linear,  we  Introduce  two  linearizing 
approximations. 

For  the  flrst  approximation,  let  t”  denote  the  mth  time  step  of  the  calculation,  when  the  solution  is 
known,  and  t*^‘  the  next  time  step,  when  the  solution  is  to  be  found.  At  r  ■  R,  we  assume  that 

|T(R,.  t*^')-T<R,.r)  |«T(R,.t*). 

Then  Equation  6c,  at  r  ■  R,  and  t  ■>  t*^‘,  is  approximated  by 

k,  dU/dr  +  [h.  +  4Fo  CT)*)  T"*'  -  h.  T^  +  Fo  (Tj  3(1*)*]  .  (7) 

This  is  a  reasonable  approximation  because  the  temperature  changes  most  slowly  at  the  outer  wall; 
fUnhermore,  the  time  increment  At  ■  t”^'  - 1"  can  be  decreased  to  ensure  the  above  inequality  condition. 
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Figured.  Thermal  Difftisivitv  vs.  U. 

The  aecond  approximation  consists  of  expressing  the  T  in  Equations  6b  and  7  as  a  linear  function 
of  U.  The  dashed  curve  in  Figure  5  shows  that  this  can  be  accomplished  with  reasonable  accuracy  by  the 
use  of  two  straight  line  segments.  Thus  T  >  biU  +  b},  where 


» 

^  ,  •  0.332223 
‘  -  0.4191495 

Since  it  is  necessary  to  choose  one  of  the  two  line  segments  a  priori  at  time  t  ■  t”^',  we  base  the  choice 
on  the  U  (inner  or  outer  wall)  at  time  t  ■  t". 

3.3  Transformed  Radial  Coordinate.  We  Introduce  a  transformation  r  ■  r(t)  for  0  ^  s  1  (Oerber 
and  Bundy  1991)  so  that  the  constant  increment  61^  will  cluster  the  nodal  points  closely  together  near  the 
inner  wall,  where  T  and  U  gradients  are  the  largest.  We  define  the  transfonnatlon  in  the  following  two 
steps: 

C  (O  ■  +  (1 -Y)  .  (0iYil,p>2), 

r  -  d;  +  R, 


if  2,159.2  .  ,  -  96.34025  if  2, 159.2 

if  U>  2,159.2  ^  --91.33484  if  U>  2, 159.2. 
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(8) 


where  D  ■  Rg  -  Ri,  and  Y  and  ^  are  chosen  constants.  We  chose  values  of  y  ^  0.092  and  ^  &  2.25,  which 
provided  a  suitable  distribution  of  nodes.  Note  that  r  ■  Rj,  R,  correspond  to  (  »  0, 1,  respectively.  The 
actual  computations  are  tiien  carried  out  in  the  t  space:  a  restatement  of  the  problem  in  1^  and  t  is 
provided  in  Appendix  A. 

4.  INPUT  DATA 

A  detailed  discussion  of  the  input  to  the  computations  is  given  in  Gerber  and  Bundy  (1991).  Briefly) 
however,  T|  is  computed  at  chosen  stations  along  the  bore  from  the  NOVA  code  (Gough  1980)  and  h, 
is  computed  from  the  Veiitay  code  (Chandra  and  Fidier  1989a,  1989b),  which  uses  T,  and  other  NOVA 
variables  to  detemaine  h,  by  the  method  of  Stratford  and  Beavers  (1961). 

Figures  7a  and  7b  show  representative  T,  and  h,  histories  at  two  stetions  on  a  120-mffl  M2S6  gun 
barrel.  It  is  seen  that  T,  snd  h,  remain  constant  until  the  base  of  the  projectile  passes  the  given  station 
at  time  t  ■  ti<  At  this  time,  these  variables  rise  suddenly  to  a  peak,  then  they  decrease  gradually,  with  h, 
decaying  significantly  faster  than  T,. 

All  the  compuutlons  reported  here  were  performed  for  a  120*mm  M2S6  tank  gun  firing  a  DM13 
round.*  The  density  of  the  gun  barrel  metal  is  taken  to  be  p  ■  7,827.0  kgAn’.  The  thermal  conductivity 
and  specific  heat  are  supplied  by  tables  that  are  plotted  in  Figures  1  and  2,  respectively.  The  value  for 
h«  ■  6.0  J/(m*  8  K)  wu  obtained  from  experiments  conducted  by  Bundy  on  a  shrouded  M256  barrel. 

5.  FINITE-DIFFERENCE  CALCULATION 

For  the  finite-difference  calculations,  the  interval  0  ^  ( si  1  (corresponding  to  R|  $  r  ^  R,)  is  divided 
by  equally  spaced  nodal  (or  grid)  points  into  NI  subintervals.  The  constant  1^  increment  is  ■  1/NI,  and 
the  location  of  the  nodes  is  given  by  ^  -  (J  -  1)  A(  (j  ■  1,  2, ...  NI  -t- 1).  Derivatives  at  node  j  are 
approximated  as  follows: 


*  Table  1  of  Qarber  and  Bundy  (1991)  daaoribai  the  shape  of  the  gun  banel. 
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[kg/(s  deg  K)] 


0U/3^j  -  (-3  Uj  +  4Uj*,  -  Uj*,)/(2  a  -  1) 


(9a) 


OU^j  -  (U,*,  -  Uj.,y(2  A5)  (j  -  2 . NI)  (9b) 

(dU/ai^j  -  (Uh  >  4Uj.,  -t-  3Uj)/(2  AO  (j  -  NI  -I*  1)  (9c) 

0*U/3V)j-(Uj-i-2Uj  +  U^,)/(A0*  0-2 . Nl).  (9d) 

If  we  let  the  time  Increment  be  At  ■  .  then,  in  the  Cnnk>Nicolion  icheme  employed  here  (6zlilk 

1968, 402)  to  obtain  the  solution  at  time  t  ■  (J  ■  2i  NI), 

(Uj-‘ -  Uj")  -  (At/2)  [(dU/at)- +  OU/9t)-‘lj  .  (10) 

By  Equation  A*4a,  Equation  10  becomes 

-  (At/2)  -  Uj*  +  (At/2)  Hj*,  (1 1) 

where 

H  ■  (a/D*)  [f,(0  3>u/ao  w  3u/ao, 
and  fi(0  and  fa(0  are  defined  in  Equations  A>3. 

The  finite-difference  appreximstlons  to  the  equations  snd  boundary  conditions  are  produced  by 
substituting  the  derivative  approximations  of  Equation  9  into  Equations  11,  A-4b,  and  A-4c,  and  then 
collecting  the  terms.  The  following  set  of  equations  for  Uj*^‘  may  then  be  written: 

NI*1 

r  A^  Uj"*'  -  d.  (n  -  1.  2 . Nl+1).  (12) 


The  coefficients  A,|  and  d,  are  given  in  Appendix  B.  The  d.*8  involve  U  values  for  the  previous  timestep 
t  X  t*.  In  most  cases,  we  have  used  NI  ■  100. 
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The  coefficients  A^^  (j  2 . NI)  are  linear  functions  of  o(U)  in  Equations  (B-3)  in  Appendix  B. 

and  are  not  predetennlned  constants.  A  successive  approximation  procedure  is  used  to  solve  Equation  12. 
First,  an  initial  estimate  is  made  for  the  a's,  e.g.,  which  are  known.  These  a's.  when 

substituted  into  the  A^'s,  make  Equation  12  a  linear  system,  which  is  solved  for  by  a  standard 

FORTRAN  routine.  The  next  a  approximation  is  ■  a  [(U,"*')'*].  Then  is  tfie  soluUon 

to  Equation  12  when  the  substituted  into  the  A^^’s;  and  (o^**‘)*'*  «  o[(U**')“J.  This 

procedure  is  repeated  until  conveigence  is  obtained;  i.e.. 


I  \<c^  (J  -  1, 2,  ...>n  +  1). 

where  Cmb,  Is  a  chosen  aibiitrarlly  small  positive  constant  Three  iteradons  have  been  required,  on  the 
average,  for  our  calculations. 

The  computer  program  contains  a  subroutlrie  prescribing  At  as  a  ftmctlon  of  t  witiiin  a  flring  cycle  (see 
Appendix  C).  The  At  is  made  sufficiently  small  eariy  in  the  cycle  to  resolve  the  highly  transient 
phenomena:  then  it  la  Increased  to  minimize  computation  time.  The  Crank*Nicolson  method  is  stable  for 
all  values  of  At,  and  there  are  no  restriodons  on  the  relative  sizes  of  At  and  A(. 

After  the  problem  is  solved  for  U,  T  is  found  by  inverting  the  Kirehlioff  transformation.  This  step 
is  accomplished  by  applying  three-point  interpolation  in  the  stored  T,  U  table. 

6.  ACCURACY  CHECK 

The  condition  of  energy  balance  can  be  employed  as  an  accuracy  check.  We  define  two  integral 
ftmetions:  (1)  Qi. «  net  energy  per  unit  length  that  has  entered  and  exited  the  barrel  since  t  ■  0,  and 
(2)  Qb  >■  change  in  barrel  heat  content  per  unit  length  since  t  ■  0. 

Qp  -  2JtR,j[‘  h,(t)[T,(T)  -  T(R,.x)]dt  -  2JtR,j['  h.[T(R,,x)  -  TJdX 
-  3.5619xlO-'’FR,j['rfr(R,.x)}^  -  Tjjdx 


J/m 


(13) 
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Q,  -  2Kpj[^‘[£“’s(f)df]  rdr 


J/m 


(14) 


A  necessary,  but  sufficient,  condition  for  accuracy  is  that  ■■  Qf  at  each  time  step. 

7.  COMPUTATIONS 

Computations  were  perfonned  to  determine  the  heating  of  a  120*mm  M256  gun  tube  firing  DM13 
rounds.  We  chose  Just  two  of  the  many  possible  firing  scenarios.  The  first  chosen  firing  sequence 
consisted  of  16  rounds  at  40-8  intervals,  followed  by  5  rounds  at  30-8  intervals.  Figure  8a  compares  the 
two  outer  wall  temperature  histories  at  z  «  4.30  m,  about  1  m  from  the  muzzle.  At  the  end  of  the  firing 
sequence  (t  •*  780  s)  there  is  a  3.8  K  discrepancy  between  the  two  computations.  This  represents  about 
2.6%  of  the  total  outer  wall  barrel  temperature  change  (total  temperature  change  is  defined  as  T[t]  -  TJ. 
A  M-soalc  plot  in  T  vs.  t  at  the  inner  wall  would  not  reveal  the  differences  between  linear  and  nonlinear 
output  Figure  8b  shows  a  greatly  magnified  picture  of  the  histories  at  one  of  the  inner  wall  temperature 
peaks,  where  the  discrepancy  maxima  occur.  The  difference  at  this  point  is  about  2S  K,  or  about  3.4% 
of  the  inner  waU  barrel  temperature  change. 

The  main  comparisons  were  made  for  our  second  firing  sequence,  namely,  the  one  shown  in  Table  1, 
which  is  almost  idendcal  tu  Scenario  #2  in  Table  1  of  Artus  and  Hasenbein  (1989).  (Our  final  burst  is 
seven  rounds;  that  of  Scenario  #2  is  six  rounds.)  This  scenario  represents  the  case  where  all  the  MIA  1 
tank-stowed  ammunition  is  fired  as  fast  as  possible,  with  two  intermediate  cooldown  periods  corresponding 
to  accessing  different  ammunition  storage  compartments.  We  have  computed  the  change  for  41  rounds, 
whereas  the  MlAl  tank  stows  40  rounds.  Also,  we  have  assumed  that  all  41  rounds  are  DM13  (kinetic 
energy  [KE]  rounds),  whereas  in  actuality  the  MlAl  stores  a  mix  of  KE  and  HEAT  (high-  explosive 
antitank)  rounds. 

The  simulations  were  performed  at  two  locadons  on  the  barrel,  4.30  m  and  0.615  m  fmm  the  breech, 
for  a  high  and  a  low  ambient  temperature:  322  K  (120"  F)  and  233  K  (-40°  F),  respectively.  These  two 
locations  correspond  to  positions  where  the  barrel  is  relatively  thin  and  thick,  respectively.  The  two 
ambient  temperatures  correspond  to  frequently  used  Army  hot  and  cold  testing  conditions.  The  initial  gun 
barrel  temperature  was  the  some  as  that  of  the  surroundings. 
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Table  1.  Firing  Sequence 


1.  17  rounds  it  7  roundiAnln 

2.  5  min  cool  down 

3.  17  rounds  at  7  iDundi/min 

4.  5  min  cool  down 

5.  7  rounds  at  7  itninds/nin 

6.  Cooldown 


Figure  9a  presents  a  comparison  of  the  two  inner  wall  heating  predictioiui  at  z  a  4.30  m  for 

m  322  K.  (For  convenience.  T  (R|)  has  been  truncated  at  600  K.)  The  largest  difference  occurs  at  the 
end  of  each  firing  bunt;  at  t  ■  1,000  s,  for  example,  the  evaluations  differ  by  about  4.8%  of  the  total 
temperature  change,  similar  to  the  results  from  our  first  firing  sequence.  Hgure  9b  shows  a  greatly 
magnified  picture  at  one  of  the  temperature  peaks.  The  relative  difference  between  maximum  values  is 
about  3.3%  of  the  total  temperature  change. 

In  Figure  10,  the  outer  wall  temperature  calculations  are  compared  for  the  two  ambient  temperatures 
at  z  ■  4,30  m.  At  t  ■  1,000  a,  shortly  after  the  end  of  the  final  bunt,  difforences  are  4.S%  and  2.5%  of 
total  temperature  change  for  T„  ■  322  K  and  T«  ■  233  K.  respectively. 

The  next  two  figures  show  heating  results  at  z  ■  0.61S  m,  close  to  the  breech,  for  the  scenario  of 
Table  1.  In  Figure  1 1,  the  linear  and  nonlinear  calculations  are  seen  to  be  almost  identical,  more  so  even 
than  at  z  ■  4,30  m.  However,  the  Inclusion  of  temperature  dependence  results  in  a  slightly  higher,  rather 
than  lower  (as  in  Figure  10),  outer  wall  temperature  prediction.  Bgure  12  depicts  the  inner  wall 
temperature  histories  under  cold  and  ..ot  firing  conditions;  the  linear  and  nonlinear  curves  are  almost 
coincident,  at  least  for  the  cooling  periods. 

Figure  13  shows  comparisons  for  the  spatial  variation  of  temperature  from  inner  to  outer  wall  for 
T«  ■  233  K.  At  z  ■  0.61S  m,  the  barrel  thickness  is  95  mm,  compared  to  23  mm  at  z  ■  4.30  m,  and  the 
temperature  is  not  able  to  equilibrate  as  thoroughly  during  the  S-mlnute  cooldown  periods  as  it  does  at 
the  thinner  cross  section.  This  behavior  is  reflected  in  the  fact  that  T  at  the  thicker  section  is  still  varying 
significantly  from  the  Itmer  to  the  outer  wall  Also,  note  that  the  differences  between  the  linear  and 
nonlinear  predictions  are  noticeably  greater  for  thinner  barrel  cross  sections  than  for  the  thicker  cross 
sections. 
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8.  DISCUSSION 


This  study  investigated  the  "error"  introduced  into  gun  barrel  heating  calculations  by  assuming 
constant  values  of  thermal  conductivity  and  specific  heat  for  the  metal  Olnear  model)  vs.  incorporating 
temperature  dependence  into  these  thermal  properties  (nonlinear  model).  The  constants  k  and  Cp  chosen 
are  within  the  range  spanned  in  the  course  of  a  typical  flring  scenario,  figures  1  and  2  show  that  the 
constant  values  differ  significantly  from  the  more  accurate  temperature*dependent  k  and  Op  at  very  high 
temperatures,  say  above  800  K.  For  tank  guns,  liowever.  limiutlons  on  the  filing  rate  and  number  of 
rounds  that  can  be  fired  greatly  restrict  the  duration  of  extremely  high  temperatures,  so  that  drastic 
differences  need  not  occur  between  predictlotui  of  the  linear  and  nonlinear  models. 

It  is  not  feasible  to  conduct  paratneter  studies  for  a  wide  variety  of  parameters  and  firing  scenarios. 
We  have  concentrated  oh  what  we  believe  to  be  the  two  cases  that  should  show  the  greatest  difference 
between  these  two  models:  the  rapid  expenditure  of  all  tank'Stowed  ammunition  under  very  hot  and  very 
cold  initial  banal  temperatures.  We  draw  several  Inferences  from  our  limited  computer  runs: 

(1)  Use  of  variable  thermal  properties  can  either  increase  or  decrease  the  banal  temperature  (e.g., 
Figure  13),  depending  on  the  case  (that  is,  on  the  totality  of  the  computer  input  parameters, 
including  T,  and  h,  histories). 

(2)  Diffeivnces  between  linear  and  nonlinear  model  temperatures  are  greater  for  thinner  barrel  cross 
sections  (e.g.,  Figures  9a,  12,  and  13).  This  may  be  due  to  tiie  more  fundamental  fact  that  thinner 
cross  sections  are  cycled  through  a  greater  temperature  range;  as  a  result,  the  effect  of 
temperature-dependent  thermal  propcnies  is  larger. 

(3)  The  differences  between  the  outputs  of  the  two  models,  constant  and  variable  theimal  properties, 
are  below  S%  for  the  cases  considered.  For  most  tank  gun  applications,  the  constant  value  (or 
linear)  model  should  be  acceptable.  In  fact,  it  might  be  preferred,  since  it  requires  roughly  one- 
third  of  the  CPU  time  of  the  nonlinear  model  in  each  computer  run. 
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9.  ADDENDUM:  HEAT  CONDUCTION  EQUATION 


In  one  dimension  the  heat  conduction  equation  can  be  understood  as  follows.  For  a  solid  material 
body,  the  heat  entering  per  unit  area,  per  unit  time  through  the  y-z  plane  at  x,  for  example,  wili  be  given 
by: 


Likewise,  the  heat  leaving  the  through  the  y*z  plane  at  x-fdx  will  be: 


^“•"■‘‘■iT'x  +  dx 


The  change  in  internal  energy  per  unit  time,  per  unit  volume,  will  be: 

-kill  +kill 

du .  ‘ifir'x  ^'x^dx 

It  ■  — - -31 - 


Expanding  the  temperature  gradient  at  x  •f  dx  about  x,  retaining  terms  up  to  first  order  in  dx,  then 
substituting  yields; 


du 

IT 


“ST" 


From  thermodynamics,  if  the  elements  boundaries  arc  fixed,  then  we  are  assured  that  u  ■  uCT),  and 

(lu(T)  .  du(T)  3T  .  ar 

~jr  "Tfiir  '•sr  ■ 

Hence,  equating  the  last  two  expressions,  yields: 


pc 


aT 

-■JT 


“"57“ 


If  these  arguments  were  repeated  in  three  dimensions,  the  results  would  lead  to  the  heat  conduction 
equation  of  Section  3.1,  viz.: 

pc^.^  -  V-{lr.7T}  . 
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APPENDIX  A: 


STATEMENT  OF  PROBLEM  IN  4,  t  VARIABLES 
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ImiNTIONALLY  LEFT  BLANK. 
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We  repeat  the  transformation  of  Eq  uation  8: 


C(K.)  Y)i^  (0<Y^ l,p>2)  r  -  DC  +  R,.  (A-1) 


Then 


dC/d{  d'C/nP  » C"  -  P(P  -  i)(i  -  r)4*-‘ 

{'(0) .  r.  x,m  c'(o  .  Y  ♦  p(i  -  Y) 

>,«C"(i)-P(P-i>(i  -Y). 


We  dellM  fi  end  t^i 


f,  -  i/(C0*. 

fa-(D/C0/(DC+R,)-C"/(C')*. 


(A-3) 


Ihe  transformed  problem  for  U  is 

3U/at  -  (a/D»)[f,(C)a*U/9C’  +  f,(C)3U/95)  ■  H(C,t)  (A-4a) 

3lJ/aC  -  [h,(t)DY/k,l  [b,U  +  b,  -  T,(t)]  at  C  -  0  (A-4b) 

[k,/(DX,))3U/3C  +  b,w,U  ■  w,  "  bjW,  atC- 1.  (A-4c) 
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where 


The  Initial  condition  is 


Wi  -  h.  +  4Fo(T*)*  ) 


'j  ■  h.T.  +  Forrj  +  3(T*)«] 


at  ^  ■  1 . 


(A-5) 


u  -  u.  -  ucTj  t  -  0, 0  s  ^  s;  1 . 


(A-6) 
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APPENDIX  B; 

COEFFICIENTS  OF  EQUATION  12 
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Intentionally  left  blank. 
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1. 


We  first  define  v,(t)  and  V}(t): 


V,  ■  2  A5  h,(t)  D  Y/ki,  v,  -  v,[bi,  -  T,(t)], 
where  bj  and  are  given  in  the  last  paragraph  of  Section  3.2.  Then. 

A„  -  -(3  +  b,v,),  A„-4,  A„  -  -Ij 

d,  -  V,  J  .  (B-l) 

At  die  outer  wall,  where  J  ■  NP  <■  NI  -f  1.  the  eoefflcienta  of  the  boundary  condition  are 
determinod  by  the  following  aequence  of  formulas: 

w,-h.  +  4FoCrNr*)* 

Wa-h.T.  +  FctT.^  +  SCTN,*/) 

Xj-Y+Pd-K) 

q,  ■  2  Ali  Wi  bj  D  X,/k, 

q,  -  2  A^i  DX,(W2  -  w,  h,)/k, 

^NP.NP-a  "  ^  •  ^NP.NP-l  ■  ^NP.NP  ■  (3  Qj)  I 

d^p  -  q,  J .  (B-2) 

We  next  define  foui  ftuictlons.  Pij.  py.  p^.  and  0  *  for  2  ^  J  ^  Nl: 

p^  -  -2  f,j/(AJi)» 
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p,j  =  [f,j/(AO*]  +  [f^/(2  A4)] 


g;  -  (ot,-/D*)  tpy  Uj.,-  +  ftu  u;  +  Uj„-]. 


where  fu  «  f|(^)  and  f^i  ■  fj  (^)  are  defined  in  Appoidlx  A.  Then 


\j.,-.AtO(,-‘py/(2D») 

(B.3a) 

V 

•3 

1 

a 

(B-3b) 

D*) 

(B.3c) 

d,  -  Uj*  +  (At/2)Gj“  . 

(B-3d) 

All  other  A^i  are  equal  to  zero. 
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APPENDIX  C: 
TIMESCALE 
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INTENTIONALLY  LEFT  BLANK. 


N 
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Here,  time  -  will  refer  to  time  within  one  firing  cycle  —  t^«  0  at  the  beginning  of  the  cycle. 
Six  constants  are  given:  t^,  t,,  t/,  ta,  At|  and  Ati .  Here,  t^  is  the  delay  time  for  the  rapid  rise  in  T, 
and  h,  fiom  initial  conditions,  and  tf  is  die  time  between  successive  firings.  The  time  increment  At  (tO 
is  given  by  the  following  fimcdon: 


At  -  t^  -  At/ 
At  -  At/ 

At  -  C,  +  Cjt^ 
At  ■  Ata 


0  i  t'  <  t^  -  At/ 
-  At,'  S  t'  <  t/ 
t/  i  t'  <  ta' 
ta'  S  t' 


where 

Ca  ■  (At/  “  At/)/(ta  “  t/)  and  C,  *  Ata“  C2ta 
(If  t'+  Ata'  >  tf,  set  At  *•  t,  -  tO  • 

A  typical  set  of  values  of  the  paiametera  would  be  the  following: 

t,'  «  0.018  s,  t/  «■  10.0  s,  At/  n  0.00025  s.  At/  ■  6.0  s. 


These  values  were  chosen  on  the  basis  of  experience  fiom  many  empirical  studies  made  with  die 
numerical  parameters. 
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Intentionally  left  blank. 
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LIST  OF  SYMBOLS 


A^i  coefflcient  in  Equation  12  for  barrel  temperature 

bi,  bi  constants  in  T  bi  U  -f  ba,  see  last  paragraph  of  Section  3.2 

c,  (T)  specific  heat  of  gun  barrel  [J/(kg  K)] 

D  ■  (R, -Ri)«  thickness  of  gun  barrel  [m.  mm] 

d.  coefScieiK  in  bant:  temperature  equations.  right*hand  side  of  Equation  12 

P  radiation  interchange  factor  between  barrel  outer  wall  and  environment,  Equation  4 

fi,  f]  given  fUnotiotu  of  Equation  A-3 

H  ftincdon  defined  in  Equation  A4 

h,  heat  transfer  coefficient  -  bote  to  gun  barrel  [J/(m*  s  K)] 

h.  heat  transfer  coefficient  •  gim  barrel  to  ambient  air  [J/(m*  s  K)] 

j  index  indicating  radial  location  of  a  nodal  point  (j  ■  1  for  inner  wail.  J  ■  NP  for  outer 

wall) 

k(T)  thermal  conductivity  of  gun  barrel  [J/(m  s  K)] 

k,  k(T,) 

NI  number  of  intervals  inR(ekr^R,andO£^£l 

NP  -NI  + 1 

Qb  increase  in  heat  content  per  unit  length  of  barrel  since  t «  0  [J/m] 

Of  net  quantity  of  heat  per  unit  length  that  has  entered  barrel  since  t  ■  0  [J/m] 

Ri.  Rg  radial  coordinates  of  inner  and  outer  walls,  respectively,  of  gun  barrel  [m,  mm] 

r  radial  coordinate  in  transverse  plaiw  [mmm]  (r «  0  at  axis  of  gun  bore) 

T  temperature  in  gun  barrel  [K] 

T,,  temperatures  in  the  bore  and  ambient  air,  respectively  [K] 

T«  -T. 

T|  prescribed  lower  integration  limit  in  definition  of  U  (Equation  5) 
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t  time  fit)m  initiation  of  first  round  [s.  ms] 

delay  time  at  given  z  for  nq)id  riie  in  T,  and  h,  [s,  mi) 

t,  time  interval  between  two  lucceiaive  rounds  [s,  ms] 

t'  time  measured  within  a  firing  cycle  [s,  ms] 

II(T)  Kitchhoff  transformation  Auiction,  defined  in  Equation  5  [K] 

Uj  ■U(ft-U-1]A5) 

u.  ■  ucrj 

z  axial  coordinate  (z  ■  0  at  breech)  [m] 

otCT)  ■  lc/(p  Cp).  thermal  diffiisivity  of  gun  barrel  [mVs] 

p,  Y  prescribed  constants  in  coordinate  transfonnation,  Equation  8 

At  computational  dine  increment  ■  (tf  -  tf)  [s] 

■  1^,  constant  step  size  in  ((0  ^  id  1) 

C  ■Y5  +  <1-TI)V 

6  azimuthal  coordinate  in  transverse  plane 

"  !)•  defined  in  Appendix  A 

(  transformed  radial  variable,  Equation  8;  Independent  variable  in  Equation  A-4a 

p  density  of  gun  banvl  metal  [kg/m*] 

o  Stefan-Boltzman  constant  ■  5.669  x  10  '*  J/(m*  s  K*) 

Supencript 

m  index  of  current  time  step,  when  solution  is  known 

m  1  index  of  next  time  step,  when  solution  is  to  be  calculated 
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